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We discuss the feasibility of a quantum nondemolition measurement (QND) of photon number based on 
cross phase modulation due to the Kerr effect in Photonic Crystal Waveguides (PCW's). In particular, we 
derive the equations for two modes propagating in PCW's and their coupling by a third order nonlinearity. The 
reduced group velocity and small cross-sectional area of the PCW lead to an enhancement of the interaction 
relative to bulk materials. We show that in principle, such experiments may be feasible with current photonic 
technologies, although they are limited by material properties. Our analysis of the propagation equations is 
sufficiently general to be applicable to the study of soliton formation, all-optical switching and can be extended 
to processes involving other orders of the nonlinearity. 



*Also at: Department of Applied Physics, Stanford University, Stanford, CA 94305; ifushman @ s tanford.edu I 
tAlso at: Department of Electrical Engineering, Stanford University, Stanford, CA 94305 



2 

In this paper we focus on the feasibihty of realizing the QND photon number measurement proposed in |ll in PCW's. QND 
measurements are important in a variety of quantum information processing techniques as well as quantum state preparation. 
Although our investigation was motivated by quantum information processing in photonic crystal (PhC) on-chip networks, a 
detector with the necessary sensitivity may prove to be valuable on its own. We show that the reduction of group velocity and 
small interaction volumes in PCWs lead to an effective enhancement of the third order nonlinearity and, theoretically, make 
experiments with high quality structures and attainable laser power levels feasible. 



We consider the case of a signal pulse from both a photon number emitter and a coherent state. A typical single photon source 

n 

is an InGaAs quantum dot (QD) coupled to a PhC cavity as in |2]. The radiative lifetime of such QD's coupled to PhC cavities 
is « 0.2 — Ins. In this experiment, a weak signal pulse is combined with a strong coherent probe in one arm of a Michaelson 
interferometer. The probe acquires a phase shift that is directly proportional to the signal photon number, and the signal pulse is 
retained for further use. The main impediment to this measurement is the small value of the nonlinearity and the relatively large 
photon absorption in semiconductor materials. 

The Kerr effect is a third order nonlinearity (x'^') and can be described by a weak intensity dependent refractive index (712/) 
as 3x^^^ = cn^n2, where c is the speed of light and n is the refractive index of the material and I is the Electric field intensity 

n 

yl. We will focus on Alluminum Gallium Arsenide (AlGaAs), which has a reasonable x^'^^ and 712 ~ 1.510 ^^^^ at a 
wavelength of 1500 nm and a high refractive index n « 3.4 ||4],|5|]. The large refractive index is attractive for the fabrication of 
PhC devices and can be combined with our current QD sources. Another material system of interest is Indium Gallium Arsenide 



Phosphide (InGaAsP), which has n2 ~ —5.9 x 10^^^^^ at 1545 nm An order of magnitude estimate of the phase due to 
the nonlinearity can be gained by expanding the index as: n « rt + n-zl- The relative dielectric constant to first order in intensity 
is then e « + 2nn2l, and the wave vector (k) in the bulk material becomes fc = fco + Afc = 7 V^- The acquired phase over a 
distance L for a signal and probe photons with wavelength A^, Xp and Ns signal photons, is then j^^j^^^PjNs, which for an area 
of A — Ifim^, length L = lOO^ni, a lifetime of — Ins and Ap_s w l.bfini gives A$ w 8 x 10^^^ x A^^. The sensitivity of 
the interferometer is given by its signal to noise ratio (SNR). The noise in interferometry with coherent states comes from the 
photon partition noise (shot noise) at the input beamsplitters of the interferometer. It can be shown that in the case of a coherent 
probe and signal in a number state, the SNR is {4(f>'^Np)~^ , where Np is the probe average photon number, and 0^ is the phase 
due to a single signal photon ItI. This means that Np > lO'^'^ probe photons are needed for the bulk experiment in order to 
overcome the shot noise when the probe is in a coherent state and the signal is in a number state. This requires a source that can 
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produce a 48k J pulse with a nanosecond width. 

PCW's offer two improvements for such a measurement. First, the group velocity of the pulse propagating in a PhC 
waveguide is reduced to u w 10^^ — 10^"^ x c 1^,0], and the local intensity increases by this factor. The effective propagation 
length of the waveguide increases as well, although the losses associated with material absorption and scattering should also 
increase with the longer effective length. This enhancement merely allows us to make smaller structures. Pulse contraction also 
means that longer probe pulses can be generated by the external pulse, and will shrink to the desired width. Secondly, the area of 
the PhC waveguide is of order (-)^. A combination of these effects leads to an overall enhancement of ( — ) « 2 x 10^, 

for a PhC area of (250nm)^ and | = 100. This reduces the energy requirement to w 10^® J in a nanosecond pulse. This 
figure does not take into account material absorption and waveguide scattering losses, which are generally the main obstacle to 
such a measurement. We show that with material absorption parameters found in literature and the probe at wavelengths above 
the half-bandgap of the semiconductors, the experiments could be attempted in a PhC device. Our analysis of the propagation 
equations is sufficiently general to be applicable to the study of soliton formation, all-optical switching and can be extended to 
processes involving other orders of the nonUnearity. 

In the theoretical proposal, the probe is sent through two arms of a Michaelson interferometer and interacts with the signal 



in one arm. The phase shift on the probe is measured on the slope of the fringe via homodyne detection 1 1]. In our case, the 
interferometer would be made in a free standing membrane of AlGaAs that is patterned by a hexagonal lattice of air holes. The 
waveguides are made by removing rows of holes. Fig. ^shows a PCW in a hexagonal lattice and the PCW dispersion for two 
modes. The group velocity is significantly reduced at the band edge {k^ = f which makes this an attractive operating point. 
Numerical precision allows us to estimate that Vg < cx 10^^. Since the Kerr effect depends on the intensity overlap, either two 
spectrally different points on the same PCW band, or on different bands can be chosen. In the first case, the intensity overlap 
is maximized, but there is a potential for a large group velocity mismatch. In the second case the mode overlap is sacrificed in 
favor of matching the group velocities. In principle, ^ w 10^^ for a Ins pulse, which means that ponts with very close j 
values can be chosen, and the proximity is limited by the ability to filter, or by the wavelength requirements for the pulse and 
probe. 

First, we derive the equations of motion for the signal and probe pulse in the interferometer. The eigenstates of the PhC 
waveguide are solutions to; 



n 
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FIG. 1: Top left: waveguide mode dispersions calculated by the 3D Finite Difference Time Domain (FDTD) method. The solid (black) line 
is the light line in the photonic crystal, above which modes are not confined by total internal reflection.Top right: group velocities of the two 
modes derived from the disperion curves by numerical differentiation {vg ~ Bottom left: modes of the PhC lattice at the fc = point. 
The mode with even vertical symmetry relative to the waveguide axis corresponds to the red dispersion curve and the odd symmetry mode is 
that of the blue dispersion curve. Bottom right: A scanning electron micrograph of a fabricated PCW in AlGaAs. 



where e(r) is the relative spatially varying waveguide dielectric constant. The solutions are Bloch modes u™(7^*)e*('^^^'^('^)*) 
where {r + az) = (r ) for the lattice with periodicity a, z is the direction of propagation along the waveguide, and m is the 
index of the band of particular symmetry. These modes satisfy the wave equation: 



V X V X K(r)e*'=^) = 



e{r)ukir)e 



ikz 



(2) 



The Bloch state is re-normalized for convenience in the last part of the paper The waveguide modes can be shown to obey the 
following orthogonality conditions : 



/.»„(.>r"j.-<'-''" = v.4. (3) 

Jn 

where the integral is taken over the whole space fi. In what follows, index m is dropped, unless it is necessary, and e = e(r). 
The waveguide modes can be rewritten to solve a different Hermitian operator: 

6 = X V X ^ (4) 

The eigenstates of this operator are {r\u, m, k) — ■yeu™(r)e*'-'^^^'^^'^-'*'' with eigenvalue , and {u, m, k\u, n, I) — 6m.nSk,i 
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by Eq.|3l The inner product denotes integration over all physical space. 

A pulse propagating in the PhC waveguide in the presence of the weak nonlinearity may be written as E = 
J dkA{k, t)\u, m, k) where A{k, t) is a time dependent coefficient of each k component. The k-space range over which 
the integrand is appreciable depends on the frequency distribution of the pulse. For pulses with ^ ~ 10^^, and group velocity 
of V = ^ sa 10^"'. So the integrand in the expression for E is dominated by a particular k about which the pulse can be 
expanded: 

expanding = cj(fco) + ^\k=ko{k - fco) + \ ^\k=ka{k - fco)^ = t^o + Vgq+ j^q^ with q = k - ko gives: 



E w Uk„if}e''-'">'-'^"'^ J + fco,t)e'['^'"''''*^"^'''^*] = (f)e*('=''"-"''*) x F{z,t) = -^|u,fco) x F{z,t) (6) 

Here F is a slowly spatially and time varying envelope of the signal or probe, which extends over many periods of the waveguide. 
In order to determine the interaction of pulses propagating in the PhC waveguide, we need to know the evolution of such an 
envelope. In the Appendix, first order perturbation theory is applied to the operator O to determine the evolution of the Fourier 
components of the envelope. To first order in the nonlinear perturbation, and with negligible group velocity dispersion, the 
evolution of two pulses (S and P) in the same waveguide, but possibly coupled to different waveguide modes (s,p) is given by 
(see Appendix): 



1 1 

S = i-tiioAlsJSl^ + 2ls,p\P\^)S - VsS' + i^-^S" (7) 

P = i^KLOp{jp,p\P\'' + 2jp^s\S\'')P - VpP' + ^\^P" (8) 

with: 

ls,p = - I d\e{T)\u,\''\up\^ (9) 

The dot in the above equation denotes differentiation in time, and the prime is a derivative in the direction of propagation (z). 
The overlap 7s, p, has dimensions of m^^, due to our re-normalization of the Bloch state to d^rjitp ~ a. The function e has 
a value of in the material and is zero in air. Here, p and s label the probe and signal modes at a particular k point, and Vs 
and Vp are the group velocities of the signal and probe pulses respectively; k is defined as cenn2 (see Appendix). The intensity 



of each pulse is enhanced by ^p-, relative to bulk, as can be shown from Poynting's theorem 



^2 (Si 



When there is no nonUnear 
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FIG. 2: Amplitude of E field for k = | ^ , | ^ and ^ 

TABLE I: values for coupling 7 for different modes of the waveguide in units of -ij, and mode volumes for each unit cell of the waveguide (in 
units of ^ . 1 and 2 refer to the first and second modes of the waveguide. 





1,1 


2,2 


1,2 




6.4 X 10"^ 


7.9 X 10"^ 


1.4 X 10"^ 




3.9 X 10"^ 


2.8 X 10"^ 


2.5 X 10"^ 



coupling, each pulse propagates with group velocity Ws,(p), and spreads according to i — g^^- The above equations are derived 
in the Appendix and can be used to investigate self focusing, soUton formation, and other effects in PCW's. In the presense of 
the nonlinearity, the pulses experience self-phase modulation due to 7s(p).s(p) and cross-phase modulation due to 7^ p terms. The 
integral for 7s^p, which gives the coupling strength, is taken over a unit cell of the waveguide, and is normalized by the length of 
the period a (see Appendix for further details). 

The shape of Us and Wp, and hence the values of the 7 terms, is not strongly k dependent within numerical error for a wide 
range of wavevectors, as determined by 3D Finite Difference Time Domain (FDTD) simulations (Figure |2}, and SfeUfc w 0. 
Thus, the coupling strength 7^ p only depends on the waveguide branch for the modes and not the particular k point. The total 
effective interaction strength is k dependent, since the group velocity determines the propagation time. The coupling strengths 
in units of a^^, and mode volumes of the unit cell of the waveguide in units of oT^, are shown in the Table |I] The mode 
volume for each waveguide mode is defined as Vi.i(2.2) ~ iA'^\{'2.)^)^%ax /a '^'^^£|"i(2) and the mode volume for the overlap 
is Vi,2 = (e|wi||ui|)„i^X\d3re|ui||wi|. 

Each equation can be transformed into a coordinate frame moving with the probe and signal respectively via x = z — Vgt and 
X — z — Vpt. The dispersion terms in0complicate the solution. We will assume that the length of the waveguide is small enough 
so that the measurement of the induced phase and the measurement of the phase on the probe is unaffected by the dispersion 
throughout the propagation. With and neglected, the solution and upper bounds on the phases on the probe after time 



P(z') = Exp[^iKojp I (7p,pl^(z')l' + 27s,p|5(2' + ^vt')\'')dt']P{Q) = p{Q)e'(^p+*s) (iq) 
Jn 

cj,p = l-KUJp[ Jp,p\P{z')\^dt' ^l-KUJpJp,p-\P{z')\^ (11) 

(j)s = nujp / ls.p\S{z' + Avt')\^dt' « Ktjp7,,p-|S'(z')P (12) 
Jo ^ 

(13) 

The phase (j)s is the phase on the probe due to the nonlinear interaction with the signal, and gives the signal photon number 
In the appendix we derive that the ideal case of a negligible group velocity mismatch and a narrow probe, gives the phase shift 
per signal photon of: 

L 1 

4>s,ideal = Cn2js,priUJsi^p (14) 

Vp VsTs 

L is the length of the PhC and is the temporal width of the signal wavepacket. 

Both the signal and probe wave undergo material absorption and scattering due to waveguide losses. PCW losses are already 
below IdB/^m and will improve with time. The material absorption consists of the linear absorption coefficient ai and the 
nonlinear coefficients, of which we will only consider the two photon absorption coefficient a2- In the case of AlGaAs at 



the half bandgap, the values of ai and a2 were found to be w .1cm ^ and « -2^^ respectively |4|]. Thus, ai limits us 
to L— sa 10cm. For a} and sub-uJ pulses, 012 results in a smaller attenuation length on the order of 50um at best. Thus, 

Vp 

experiments with the signal and probe at the half-gap are not feasible. In order to circumvent pump depletion due to two photon 

J 1 

absorption, a pump at even longer wavelengths above 1550 nm should be used II IB . In that case, 02 is close to zero, and we will 
assume that the 100 /im PCW length is the limit. In this case, the pump will propagate in the lower branch of the waveguide, 
while the signal should couple to the upper branch. For example for a signal at 1550 nm in the upper waveguide branch, a pump 
at 1620 nm should be used in order to have both beams velocity matched at the vr/a point. We briefly mention that the GalnAsP 
material system has ai « IcmT^, and n2 w 5 x 10^^^^^ and most likely similar two-photon absorption, which means that 
both materials are suitable candidates for an experiment. While the nonlinearity is enhanced closer to the band-edge of the 
semiconductor bandgap, the absorption increases accordingly and reduced the interaction length. 



The phase due to a single photon in signal S and the energy required for an SNR of 1 for number state detection in AlGaAs, 
are plotted in Fig. |3] We plot both the ideal case, in which two-photon absorption is negligible, and the reality in which the 
pump is at 1620nm. 
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FIG. 3: Phase shift due to a single signal photon with a lifetime of 200 ps, after propagation through a 100 AlGaAs waveguide with a 
narrow probe and no group velocity mismatch (A), and the energy required for an external pulse to obtain a SNR of 1 (B). In (A) and (B) it 
is assumed that the signal and probe are at 1500 nm and two-photon absorption is not present. In (C) we plot the phase for the case of the 
signal photon in waveguide 1 at 1550 nm and probe at 1620 nm in waveguide 0. The required probe energy for this scheme is shown in (D). In 
all plots, the blue and red curves correspond to both the signal and the probe in waveguide modes or 1. The black curve corresponds to the 
probe and signal in different waveguide modes 

There are two sources of noise in this experiment in the case of an ideal detector. One is the phase noise due to intrinsic noise 
of the signal beam, and the other is the interferometer noise due to the uncertainty of the probe photon number Following O], 
it can be shown that in the case of a coherent signal state with mean photon number (n^) = Ng and coherent probe with mean 

s^intrinsicl 02^2 



photon number {fip) = Ns the uncertainty in the detected signal is (An^ ^^3^^^,^^) = (An^ j^j^^^^j^) + ^^jfl- There are two 



cases of interest: the signal in coherent and number states. For the coherent state, 

1 1 



^(t>lNp (15) 



i^^l, observed) 

In the case of the signal in the number state, the intrinsic noise of the signal disappears, 

{^''^l, observed) — . ,2 AT ^^^^ 
WsNp 

When the probe photon number is reasonably large, we can relax the requirement on Np. If the tolerated error for coherent state 
detection is _E = PNg, then the condition is Np = (40^/3iVs)^^, and j3 < 1. Thus, detection of 1000 signal photons with an 
error of 100 (/3 = 0.1), would require 50 — lOOnJ. For smaller signal photon numbers, the level of tolerated error decreases, 
and the requirement is more stringent than number state detection, since > 1. 



In conclusion, we have derived the equations of motion for a probe and signal wave interacting via the third order nonlinearity 
in a photonic crystal waveguide. Within the slowly varying envelope approximation, the equations yield intuitive results, and are 
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essentially identical to the equations of propagation for pulses in nonlinear fibers and materials, if the plane waves used in the 
mode expansion of the electromagnetic fields are replaced by Bloch waves. However, the user of PCW leads to the necessary 
enhancement of the pulse intensities due to the small mode volume and reduced group velocity of the pulses. We have shown 
that for the case of a very long wavelength probe pulse, that does not suffer from two-photon absorption in the AlGaAs material 
system, the energy requirement on the probe wave is within attainable values (w in sub ns pulses). Since the sources of such 
pulses are external to the PCW, the generated probe pulse can be broader than the pulse desired in the PCW, due to contraction 
by the group velocity. Our derivation has assumed that coupling into the waveguides and the beamsplitter implementation in 
a PCW are perfect, and the scattering loss of the waveguide can be neglected. This is of course a gross generalization. High 
couphng efficiencies and low loss propagation over lO's of microns have been shown, and will only improve in time. We are in 
general most strongly limited by two photon absorption in this proposal. Other material systems such as GalnAsP Q], which 
also exhibit higher n2 values that AlGaAs, could also be considered for this implementation, but ultimately atomic resonance 
systems and systems with a high phase shift and low loss are necessary | ^ 13 1. In principle, an on-chip QND photon number 
detector could be a component of a photonic crystal based quantum circuit, or can serve as a sensitive intensity detector and 
switch. The derivation presented here, can be easily extended to other types of intensity and field dependent nonlinearities, and 
can be used to analyze other nonlinear optical effects in PCW's, as well as soliton formation and propagation. 
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APPENDIX A: DERIVATION OF THE PROPAGATION EQUATION 



The probe envelope evolves according to (with q = k — ko): 

= ^;,^(/,^i)ei[(fc-feo)^-(a;(fe)-a;o)t] 



dP d 
'di 

« 11 dqAiq + ko,t)ei''^'-'''''^-'''''^'- 

1 dv. 



J ^^^ dA{q + k„t) _ ^^^^^ ^ ^^q^)A{q + ko,t))e 



(Al) 
(A2) 
(A3) 



= / ''^ '^^^'ot '°'^^ e4^(---*)-^^-^*] - (.,| - 4^^) J dqAiq + ko,t))ei^^-^^^-'^^'^'] (A4) 



in the case of nonlinearity, the wave equation with c 



Moeo is: 



-VxVxE= — ' + fj,o- 



9*2 



And we can rewrite this equation in terms of the previously introduced Hermitian field operator as: 



-^V X V X -^^^eE = -O^eE = + /xo^^^ 



9*2 



(A5) 



(A6) 



(A7) 



Here, in general, the polarizability P is given by Pi = eo Ejvt,; Em,n,p x\j,k,iE{'^p)iE{uJs)jE{u/p)k, where xlj,fc,i ^'^^ 
components of the third order nonUnearity tensor. In our case, we consider the system to be isotropic, the response instantaneous, 
and only two frequencies (Wg, tOp) to be present. To find the time evolution of the coefficients, we use the wave equation. In the 
presence of the nonlinearity we expand e ~ e + (5 and O as: 



6 + AO = -^V X V X ^ - - 
Ve ve 2 



S 1 ^ ^ 1 1 ^ ^ 1 ^ 
— X V X ^ + X V X —=- 
e ve ve V^ ^ 







+ 


W] 







o + A6«d-i[^o + 6^] = o-i{J,6} 

Let {u'\ = {u, k', n'\, \u) = \u, k, n), A = A{k, i). Then we have: 



(w'|(6 + a6) J dkA\u) = -{u'\^^ JdkA\u) 

j dkA{{u'\d\u) + {u'\Ad)\u) = -{u'\^ j dk{A- cu^A- 2iioA)\u) 

J dkA{u'\Ad\u) = ^ y dkA{u'\u) 

-i / dkA{u'\-d + d-\u) = ^A{k') 
2 J e e c 



.1 I dkA{ur-\u){^ + '^)^^A{k') 
-'^ [ dkA{k){u'\-\u) = ^i(fc') 

C J 6 C 



(A8) 
(A9) 

(AlO) 
(All) 
(A12) 
(A13) 
(A14) 
(A15) 
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Above, we neglect second derivatives of the envelope and combine the two frequency terms as w ^ + w ^uP^ + Vs{k — k') w 
2w^, because the (fc — k') term will lead to the derivative of the slowly varying envelope multiplied by the nonlinearity and is 
very small. Since the frequency bandwidth of the envelope is small, the envelope is slowly varying in time, and the second order 
time variation in the coefficients A{k) is neglected. Furthermore, we have also assumed that the inner product {u'\u) is roughly 
unchanged by the nonlinearity - it remains a delta function. The perturbation Sg.p contains both the real and imaginary parts of 
the third order susceptibility. The real part is responsible for the cross phase modulation, while the nonlinear term gives the two 
photon absorption of the signal and probe. The figure of merit for the feasibility of the experiments is the phase-shift gained per 
loss length (-j point). The full perturbation can be written as: 

5s,p = mi + 3eo(xl'^ + ixf + 2|^K,.)n (A16) 

Where x'^' is the third order polarizability (which is assumed to have only one value and to be infinitely fast), and E{ijJs,p) — 
{5, P}ufc^ ^gi(fcs,p2-i^a,p*) jg j-jjg electric field of the two modes. The linear loss ai is the imaginary part of the dielectric 
constant. 

(3) 

The value of Xr can be determined from the experimentally observed bulk intensity dependent refractive index defined via 
h — n + n2l, where I is the average field intensity. For a bulk material 1 is given by / = = iceon|i?p. Thus, 

3^(3) _ cn^n2- For AlGaAs at the wavelength of 1.5 /im, n2 ~ 1.5 x 10~^^^^; furthermore, the index is very similar for 
TE and TM polarization in AlGaAs slab waveguides j4l. Since x*^'^^ measures the response of the local charge distribution to 
the local electric field, the coefficient itself is a material property and is not modified in the PCW, except possibly due to surface 
effects (e.g. reduced response or artificially added birefringence). Thus, we can derive the value of the coefficient from bulk 
experiments and combine it with the modified electromagnetic fields to get the resulting effect in the PCW. The x^'^^ coupling 
term only exists in the material, and we can replace the term with a dielectric which is equal to the spatially patterned index 
of AlGaAs in the PCW and is zero in the air. We set Sx'-'^-' = cn2e{r)- And we define k — ceof^2, so that the perturbation due 
to real part of the nonlinearity becomes ds,p = Ke{r'}{\E{uJs,p)\'^ + '^\E{ujp^s)\'^) ■ The loss terms are similarly determined from 
a fit to atotai — ai + ct2l- The linear loss ai results in an exponential decay of the signal with a characteristic length (ai)^^. 
The nonlinear loss gives a characteristic length of (0:2/)^^. We will drop the losses for now, in order to derive the equation of 
motion for the pulses, and will assume that the Bloch components of the eigenstate \u) and \u') belong to the same waveguide 
branch n = n' . Each of Un,k, Un,k' is then roughly given by some central k component that is modulated by an envelope 
Un,ke^''^ ~ y^^^gi(feoz-c^ot)g«[(fe-feo)z-(t^(fe)-t^o)t]^ \Yg jjQyy insert the exact form for and \u) into the above equation, and 
only look at the cross phase modulation component on the probe due to the signal. We take n to be branch of the pump mode p, 
and m to be that of the signal s, and only treat the perturbation due to the signal explicitly. Also, we will assume for simplicity 
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that the signal group velocity at k and k', as well as the dispersion so that we can combine them in the expansion of the 
Bloch state. Eq. lAlSl then imphes: 



i(fc') = liupj dk" j dV'A(fc")«;e>pnu,|2|5(z')|2e*P"-'=')-'-(-('=")— (*'))*] (A17) 

«zWpK j dk"A{k") j dz'|5(z')|'e*'^'=""'''^"'"''"'''"^""^'='»*^^ j d\~e\us\''\up\'' (A18) 
KiujpK-fs,p I dk"A{k") f rf/|5(/)|2e»[(fe"-fc')2'-(-(fc")--(fe'))t] (A19) 



This is now inserted back into the evolution equation for the envelope IA5I and we re-substitute fc — fcp ~ q and keep the 
frequency term in the form uj{k) — uj{ko) for convenience. 

J dkA{k, i)e'('=-feo)2-(^(fe)-c.o)t) ^ (A20) 

= iujpK-fs^pJ dk J dk"A{k") J d/|5'(2')|2e^P"-fc)^'-(-'(fc")-c.(fc))t]g»[(fe-feo)^-(a'(fe)-'^o)t] (A21) 

^iLOpKjs^pJ dz' J dfc"A(fc")|S'(z')pe'l('="^'-''°")-('^(''")-'^°)*l) J dfce*'=(^^^') (A22) 

= iLUpK-f,^pJ dz'\S{z')\^e''">^'-'-^'^S{z- z') J dfc"A(fc")e'[''="~''''^^'"("'(''"'~'^°^*l^ (A23) 

= iLUpKjs,pJ dz'\S{z')\^e''"'^'-'"^6{z - z')P{z') (A24) 

= iLjpKj,JS{z)\^P{z) (A25) 

The ^/e terms cancel the denominator of the perturbation. The term js,p contains an effective area integral 7s p = 
a /a '^'^'"^I'^sPl'^pP ~ J dxdye\us\'^\up\'^ . Since the Bloch states are periodic, their integral in each unit cell (js.p) is the 
same, and we simply weigh it by the average value of the slowly varying envelope in that cell to find the integral over the whole 
volume. A is the unit cell volume. S and P are only functions of time and the propagation coordinate (z here), and are uniform 
in the transverse (x,y) plane. The term k contains the strength of the nonUnearity. We now re-normalize the Bloch state and the 
field: 



N.hio, = J J J dVeoe(r)|5nu,|2 (A26) 

« J dzeo\S\^^ J^d^re{r)\u,f {All) 

/....ISP ,A28) 

l = -f d^fe{r)\u,\^ (A29) 



a Ja 
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The above normalization means that [7] = and [\S\'^] = Volt^. Thus the term K,'--fuj\S\'^ has units of s~^, as desired. We 
now insert the form for A into I A5 1 

The refractive index due to cross-phase modulation is twice that of self phase modulation jj]. The evolution of the slowly 
varying envelope is given by: 



1 1 

S = z-«;w,(7,,,|5p + 2-fs,p\P\'')S ~ VsS' + i^-^S" (A30) 

1 1 dv 

P = i-Ka;p(7p,p|P|2 + 27p,,|^|2)P - v^P' + ^i^^P" (A31) 

We will now show the qualitative behavior of the two pulses, assuming a weak interaction. Take P = p{z, i)e*"^(^'*) and 
S — r]{z, t)e^^^^''^\ Inserting P into lA3ll the real and imaginary parts satisfy: 
Real: 

dv , , ^ „ 
P+ivp + -g§^)P P (A32) 

Imaginary: 

+ Vp<j)'p = -KUpi-fp.pp'' + 2jsW)P + -^ip" - i.4^fp) (A33) 

If we assume that in lA32l the second derivative term vanishes, then we can see that the envelope moves along a characteristic 
given by Wp + ^^(/> « Vp{kQ + 4) ), which, in the flatterregions of the dispersion curve is very close to ^^(fco). If ^^(^) << 1> 
then IA3 3l simplifies to give an equation for the phase along a characteristic given by fp(fco + 

<i> + vp(j)' + = -KWp(7p,pp2 + 27,,p772) (A34) 

Since we will not generally be able to generate a pulse that is a solution to the nonlinear system, we will assume that the pulses 
are Gaussian and drop the dispersion terms for simplicity of the analysis. We are essentially assuming that the envelopes are 
very slowly varying and that the interaction time and dispersion are not strong enough to affect the phase measurement, which is 
dominated by the group velocity term rather than the group velocity dispersion term. Thus, we will rewrite Eq. IA31l in a frame 
moving with the group velocity in terms of z' = z — Vpt and t' ~ t, and neglect the terms with and Since the two 
pulses may have different group velocities, we have Av — {vp — Vs)- If the group velocity dispersion terms are neglected, each 
envelope is only a function of z' . 

« '-Ku:p{^p,p\P{z')\' + 27,,p|5(/ + Avt)\')P{z') (A35) 
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The solution, and upper bounds on the phases after time t = ^ are: 

P{z') = P{0)Exp['-KUJp [\%,p\P{z')\^ + 2js,p\S{z' + Avt')\^)dt'] = P{0)e^(^-+^s) (A36) 
^ Jo 

cI>p^Iko.-p [ 7p,p|P(^')l'«^*'~ ^«^'p7p,p — l-P(^')l' (A37) 

Z Jq Z Vp 

ft T 

4>s = i^p / lsAS{z' + Avt')\^dt' « Ka;p7,,p — (A38) 
Jo '^p 

If the second arm of the interferometer is adjusted for a phase shift of ■jt/2, the difference in the intensity signal on the 
two detectors gives the phase, and thus an estimate of the photon number. The total integrated signal energy is Idet = 



J J J d'^reoe\P{z')\''\un\^sin{^s)- Starting from iV.^iw, = J^^dz'eo\S{z')\^ = J^^dz'eo\S\^s{z') with \S\' 



(similarly, |P(^')P = \P?P{z') and jPp = the integral for I^et is re-ordered again: 

Idet^ J J J d''reo\P{z')fMz')e\unf (A39) 

« / dz'eo\P{z')\''Mz')- I d\e\up\^ (A40) 
J o 



/ 



dz'eo\P{z')\''<l>s{z') (A41) 

L Ns hid s flUJ-o 



/ dz'p{z')s{z') (A42) 

'eff,P J 



Vp eoLeff,S Le 

In the case when P is much narrower than the signal pulse S, we can view it as a quasi delta function p{z') « ie/ f,p^{x — xq). 
The effective length can be approximated by the spatial width of the pulse, which is TsVs, the product of the temporal width and 
the group velocity of S. In this approximation the detected intensity is : 

L NgfiLOg 

Idet ~ KiVplp,a NpTlUJp (A43) 

Vp eoTgVs 

= cn2h^colcJslp,s — (A44) 

^ Vp TsVs 

And the phase shift per photon of S is jVpi^Jfto-p ~ '^''^'^^p^s^s^p^:;;^- In order to make a comparison to our initial plane 
wave argument, we can identify jp^s as the effective inverse area, as the effective pulse bandwidth in the waveguide, and 
— as the enhanced interaction length (time). 
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